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In this work we discuss the formation of zero energy vortex and chiral edge modes in a fermionic represen- 
tation of the Kitaev honeycomb model. We introduce the representation and show how the associated Jordan- 
Wigner procedure naturally defines the so called branch cuts that connect the topological vortex excitations. 
Using this notion of the branch cuts we show how to, in the non-Abelian phase of the model, describe the Ma- 
jorana zero mode structure associated with vortex excitations. Furthermore we show how, by intersecting the 
edges between Abelian and non-Abelian domains , the branch cuts dictate the character of the chiral edge modes. 
In particular we will see in what situations the exact zero energy Majorana edge modes exist. On a cylinder, 
and for the particular instances where the Abelian phase of the model is the full vacuum, we have been able to 
exactly solve for the systems edge energy eigensolutions and derive a recursive formula that exactly describes 
the edge mode structure. Penetration depth is also calculated and shown to be dependent on the momentum of 
the edge mode. These solutions also describe the overall character of the fully open non-Abelian domain and 
are excellent approximations at moderate distances from the corners. 

PACS numbers: 05.30.Pr, 75.10.Jm, 03.65.Vf 



I. INTRODUCTION 

Models that display spinless p-wave pairing are known to 
exist in both Abelian and non-Abelian topological phases. 
The systems are BdG (Bogoliubov de-Gennes) type topologi- 
cal insulators \ and therefore support gapless chiral modes at 
the edges between Abelian and non-Abelian domains. When 
these edge modes have zero energy they are known to be 
Majorana fermions. In addition to this the bulk of a non- 
Abelian phase is capable of supporting Majorana zero modes 
which are localized, gapped, and give rise to non-Abelian 
statistics^ ^. 

The understanding of these properties has been greatly en- 
hanced through the use of exactly or nearly solvable spin mod- 
els. Arguably the most important for the spinless p-wave sys- 
tem is the Kitaev Honeycomb systenP. The Abelian phase of 
model can be analyzed using perturbation theory'^'^and is re- 
duced to the so called 'Toric Code' system in this limit^^. The 
main advantage of this s ystem however is that it can be under- 
stood as either MajorancP^^^ or Dirac fermions ^^ ^^ hoping 
in a Z2 gauge field. In the Dirac fermion picture, obtained 
using Jordan- Wigner type fermionization procedures, the spin 
Hamiltonian in each gauge sector reduces exactly to a mean- 
field type p-wave systenrP^ but where the fermionic vacuum 
is exactly that of the 'Toric Code'^^ In these paper we will 
discuss the structure of vortex and edge zero modes for the 
honeycomb system using this later representation. 

The overall aim of this paper is to present an alternative ex- 
planation for the zero and low-energy chiral modes that exist 
in this system. Our perspective is complementary to previous 
work on the honeycomb model edge states (see for example 
Ref 7 Appendix B and Ref. J/7j, the continuum p-wave anal- 
ysis of Refs. |2]and|6l and the bosonic condensation theory 
presented in Refs. '2?and[23l 

The first half of the paper describes how the notion of a 
branch cuts arise naturally from the 2-D Jordan- Wigner pro- 
cedure. This is in contrast with mean field p-wave analysis 




= 1 \Xq = -1 Vortices at the end of brach cuts 

FIG. 1 . Vortices always appear at the end of a branch cut. Figures 
(a) and (b) are real vortex configurations. Note that the eigenval- 
ues of the homologically non-trivial symmetries dictate which vor- 
tices are connected to each other. In the absence of any vortices the 
X and y anti-periodic homological conditions are encoded as lines 
X(^N:^,y) = — 1 and y(Ar^,x) = — 1 Vx respectively. With the 
conventions used irP'the term Yj^ATy ,0) = — 1 dictates which vortices 
are connected by the branch cuts. Figures (c) and (d) are "simulated 
" vortex configurations obtained by varying couplings Jx and Jy and 
hi. 



where the branch cuts are an afterthought to ensure that the 
modes are single- valued. We will see that, as expected, these 
branch cuts connect the topological defects (vortices) of the 
system. However, through out this story we will attempt to 
emphasize that it is the branch cuts that are the fundamental 
objects. For example, it is the branch cuts, and not the vor- 
tices, that dictate the fermionic behavior of the system. This 
perspective also holds on the boundaries between Abelian and 
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non-abelian domains. For example we will see that it is the 
number of branch cuts through those edges that dictate the 
character of the modes found there. In the second half of the 
paper we will analyse the zero energy bulk modes and the 
zero energy and chiral edge modes found in the model. Our 
analysis of edge modes is valid for both cylindrical and fully 
open boundary conditions but is based on the consistency rela- 
tions between homologically trivial excitations (vortices) and 
the homologically non-trivial excitations on a torus^^ . We first 
introduce the cylindrical system and describe the general char- 
acter of the modes found in this case. The general conclusion 
is that exact zero modes only form on edges that are inter- 
sected by an even number of branch cuts. In addition to this 
we see for the hard boundary condition (i.e. where the Abelian 
domain is exactly the vacuum), that there are exact solutions 
for the BdG equations. We use these solutions to examine 
the mode penetration depth as a function of the Hamiltonian 
parameters and the mode momenta along the edge. 

We finally extend the general analysis to fully open rectan- 
gular boundary conditions we see that exact zero modes only 
form in this case when there is an odd number of branch cuts 
through the domain. The reason for the difference is an extra 
phase factor that is contributed at the comers of the system. At 
moderate distances from the corners however the exact solu- 
tions for the cylindrical hard boundary system are an excellent 
approximation for the open system eigenmodes. These results 
are general agreement with Ref . 7 and Ref . 2x 



II. FERMIONIC FORMULATION 

It was shown iiPthat each vortex sector of the honeycomb 
lattice model can be written as 

q I 

where in terms of fermions we can write 

Ho = Jx'^Xq{cl-Cq){cl^^Cq^) 
Q 

+ J,^(24c,-/), (1) 
q 

where we have used the shorthand q ^= q-\-nx, q t= Q-^^y 
and q q -\- Uy -\- rix. In the plane, Yq = I for all q and 
is defined as 

y-i 
y'=o 

The terms P^^^ are explicit T- symmetry breaking terms, the 
fermionic form of which was also derived in For simplic- 
ity in this work we will retain only terms P*^^^ , P*^^^ , P*^^^ and 



P'^^\ These terms are sufficient to generate the required non- 
abelian phase and lead to more symmetrical solutions. Explic- 
itly these terms arJ^ 

= -iKX,(ct - c,)(ct^ - c,^) (3) 

P(2) = -inXg^icl^ + Cq^)(4^ + Cg^) (4) 

P(3) = +inY^icl - Cg)(4 - Cg^) (5) 

P(4) = +iKYg^{cl^ + c,^)(4^ + c,^) (6) 

The Jordan- Wigner convention used to define the fermions 
is directly responsible for how vorticity is encoded in the 
fermionic system. For the string convention chosen irPI' the 
vorticity is encoded in the fermionic Hamiltonian through the 
condition Q. On a torus there are additional homologically 
non-trivial degrees of freedom which also need to be deter- 
mined consistently with the condition (|2]). These homolog- 
ically non-trivial are encoded the Xq and Yq values at the 
boundary of the system^ ^ Recently we have extended this 
Jordan-Wigner method to deal with the Yao-Kivelson 3-12 lat- 
tice variant of the modeP 

The consistency relations provided in Ref. [211 have an in- 
teresting pictorial representation which leads us naturally to 
the concept of branch cuts and a less restrictive understand- 
ing of vorticity. For any vortex arrangement we see that there 
are lines of Xq = — 1 and Yq = —1 which together connect 
vortices in pairs. In Figure [T] we have provided a number of 
examples. 

On an open plane we no longer have these homologically 
non-trivial symmetries but neither do we have the condition 
that vortices are created in pairs: Yl^ Wq = 1. In this case 
valid vortex sectors can be encoded using the following guide- 
lines. 

• The vortex free sector (Wq = IVgr) is encoded as Xq = 
IWq. 

• A single isolated vortex at position q is encoded with 
Xq = 1 everywhere except for a single line of Xx^y = 
— 1 starting aty -\-ny and extending to infinity. 

• When two vortices occur at different x-positions there 
are two unique strands of Xq = — 1 connecting them 
both to infinity. 

• If two vortices occur at different ^/-positions but with 
the same x a line of Xq = — 1 connects them together. 

One can 'simulate' the change of vortex sectors by altering 
the coupling constants (the Jx and Jy) on unique links^^. Thus 
by changing the sign of Jx at q one effectively changes the 
gauge encoding Xq. Strictly speaking this does not change the 
vortex sector of the Hamiltonian however. With our fermion- 
ization convention, and on a plane, there is no vortex sector 
which would correspond to the change Jy ^ — Jy at q. 

From now on we will take J = Jx = Jy, dropping the sub- 
script and take the viewpoint used in^^ where, by changing 
the coupling strengths, we can simulate changing the vortex 
configurations. In what follows however, and only for conve- 
nience, we will generally continue to regard the J and k, terms 
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as constant across the lattice and allow vorticity to be encoded 
in the X and Y terms. With this perspective it is easier to ap- 
preciate that truly meaningful objects in this story are not the 
vortices themselves but the connected strings of — I's defined 
on the Xq and Yq matrices. Indeed as we have already shown 
these strings take on the role of branch cuts in our fermionic 
Hamiltonian and will see later that it is their ends that give 
rise to localized zero modes. From this perspective we can 
say that zero modes are only associated with vortices because 
a branch cut always happens to end there. 

In addition to the vortex zero-modes we will also see in 
what follows that it is the branch cuts that are directly re- 
sponsible for the appearance of the single extended zero mode 
that occurs at the interface between abelian and non-abelian 
phases when an odd number of (ordinary localized) zero- 
modes are in the non-abelian bulk. The parameter J dictates 
which phase we are in. For J < 7^/2 we are in the abelian 
phase and for J > Jz/2 we are in the non- Abelian phase if 
K, ^ 0. In what follows we will specify the J and k, values in 
the Abelian domains as J a and k,a respectively. 



III. BULK MAJORANA FERMION ZERO MODES 



have 2N zero energy {E = 0) fermionic modes of which N 
must be identified as a^'s and N as a's^^. It is rather remark- 
able that one can always choose a superposition of the 2N 
and a zero-modes such that the resulting modes are fully 
localized around the vortex excitations. 



N 



= [c|,...,cj^, Ci,...,Cm] 



(12) 
(13) 



It is interesting to note that this localization condition also 
enforces the condition that Un 



gzO^^*^ However, if one 



wishes to call this a Majorana mode = 7] it is necessary to 
multiply the states (li, v)^ by the overall phase e~*^^/^ such 
that u 



Vqj. It was pointed out by Stone and ChunJ^ 



that the Majorana condition therefore fixes the global phase 
of the states up to an overall sign ±1. Understanding how 
and when this overall sign changes is crucial to understanding 
how non-abelian statistics arise in this degenerate subspace. 
In Figurel2]we show the \uq \ and \vq \ position space structure 
for some different values of J and tz. 



In this section we will briefly discuss the bulk Majorana 
modes found at the end of the branch-cuts. We will not how- 
ever discuss the detailed structure of the bulk modes other 
than to present some numerical calucations. In later sections 
however we will demonstrate how the structure can be seen 
as a limiting case of edge modes found between domains of 
Abelian and non- Abelian topological phase. 

We begin by presenting the Bogoliubov-De Gennes formal- 
ism. The full position space Hamiltonian can be written in the 
form 



qq' 



Cq Cq 



^qq 
At 
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This system can be diagonalized by 
Bogoliubov-De Gennes eigenvalue problem 
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IV. UNPAIRED MAJORANA MODES AND EDGE STATES 

The non- Abelian phase of Kitaev models are Topological 
insulators of the BdG class, see for example^ Roughly speak- 
ing this means that we have a bulk energy spectrum which is 
'insulating' (does not cross the Fermi energy at E = 0) and 
an edge spectrum which is 'conducting' (does cross the Fermi 
energy at £^ = 0). For a careful choice of edge conditions it is 
possible to analytically treat the 'conducting' edge modes. 

In order to determine the structure of the modes let us con- 
sider an element of an arbitrary eigenstate a}^ of the BdG 
Hamiltonian. Each value Ux^y is connected to Ux±i^y and 
Ux,y±i through the non-zero elements of the ^ matrix and to 
'^x±i,y and Vx^y±i through the non-zero elements of the A 
matrix. It is quite difficult to say anything generic about the 
form that an eigenvector should have. One feature is universal 
(8) however. We see that if the elements around the point Ux^y are 
almost zero then Ux^y should also be almost zero. It is true 



(7) 



where the non-zero entries of the diagonal matrix Enm = 
En^nm ^rc the the quasi-particle excitation energies. The 
Bogoliubov-Valatin quasi-particle excitations are 



aj,...,aj^, ai,...,aM 



[ cj, cj^, ci, Cm 



U 
V 



which after inversion and substitution into ^ give 

M , 



(9) 
(10) 

(11) 



In spinless p-wave systems it is guarunteed by an index the- 
orem that in the the case of the 2N well separated vortices we 



regardless of the values we give our coefficients in our Hamil- 
tonian and it is this rule that determines the vast majority of 
the zero-mode structure (or lack of it). 

In the absence of branch cuts, there is a simple condition 
that the nine interconnected elements must obey if they are to 
be eigenstates of the system: 

(2 J, - E)ux^y + (14) 

J{Ux-\-l,y + Ux-l,y + Ux,y-\-l + Ux,y-l) + 
(J - 2ii^)vx+i,y + {-J + 2if^)Vx-l,y + 

(J + 2itv)Vx^y-\-l + { — J — 2iHi)Vx^y-l = 0. 

For edge states on a cylinder we make the reasonable as- 
sumption is that, in the direction of edge, our modes are plane 
waves (momentum eigenstates). For example along the lower 
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(a) 



(b) 



FIG. 2. (Color online) The position space structure \uq 
= IJ = 0.8 and k = 0.2 



\vq\ of vortex Majorana zero-modes for (a) = 1, J = 1 and = 0.5 and (b) 



edge of a cylindrical non-Abelian domain we have BdG exci- 
tations of the form 



AT ^ e±*^--(ii(^ - yo)4 + v{y - yo)c^) (15) 



This state corresponds to a superposition of left (right) moving 
particles and right (left) moving holes. On a cylinder the al- 
lowed values of kx are 2niT /Nx when there is an even number 
of branch cuts through the edge and 2{n ^ 1/2)11 /Nx when 
the number is odd. The basic reasoning is this. A branch cut 
is accommodated in ( p3] ) by a change in signs of the elements 
J and n acting on some (not all) of the values u and v. To 
keep the energy low then the phase of the mode a\ should 
abruptly change sign to counteract the sudden sign change in 
the fermionic Hamiltonian. 

On a cylinder this has interesting consequences. Let us start 
from the toroidal case and open up the ^-boundary above and 
below the = line. We now have two edges which are 
some distance apart. Translation invariance remains in the x- 
direction but is broken in the ^/-direction. Recall now that 
the anti-periodic x-boundary condition is encoded as a single 
line of Xq = —1. Thus in the periodic vortex free sector 
we therefore have chiral edge states with kx = ±2n7r/Nx. 
This includes two edge zero-modes, one on each edge. In the 
anti-periodic vortex free case we have no zero modes. This is 
because we have a single branch cut intersecting both edges 
and thus kx = ±2{n + l/2)7r/Nx. 

If a single vortex exists inside the cylinder there must be a 
branch cut connecting it to either infinity or some other vor- 
tex outside the cylinder. If we were originally in the periodic 
system then the introduction of a branch cut through one wall 
would destroy periodicity on this edge and we could not have 
Majorana zero modes. The other edge however would remain 
unaffected. In the opposite sense if we were originally in the 
anti-periodic sector then the introduction of a vortex would 
restore periodicity to one of the edges and thus allow values 
of kx = ±2n7r/Nx to propagate along this wall. 

We can extend this reasoning to deal with fully open bound- 
aries (non-Abelian domains within Abelian domains and vice 
versa). However it is useful to first solve the system exactly on 
a hard interface Ja = where the Abelian side of the edge is 



the full vacuum. In this scenario numerical calculation shows 
that all low-energy modes satisfy Uq = e^^Vq. Thus for modes 
along the lower edge aty = yo we have 



iO/2 J 



,+i6'/2 



(16) 



Note that under the conditions kx = and Im{f) = this 
ansatz is already a Majorana fermion. If one now substitutes 



this expression into ( 15 ) we observe that 



E{J,n,kx) 



sin kr 



(17) 



and that, along the bottom edge, 6 = tan~^(2/^/J). Further- 
more one sees that the function / follows from the recursive 
relation 



f{yn+2) 



1 



[dif{yn+i)^d2f{yn)] (18) 



where 



J2 _ 4^2 

di — 2Jz + 2J cos(kx) — i2 sm(kx) 



d2 = V J2 + 4/^2 ^ J 

Interestingly the structure of the mode depends on the pa- 
rameter Jz but the associated energy does not. However this 
feature is present for the (J a = 0) hard boundary condi- 
tion only. Indeed numerical calculation shows that even the 
sm{kx) dependence is not exact once the hard boundary con- 
dition is relaxed (Ja ^ 0). 

The mode penetration depth can be calculated easily from 
the recursive relationship ([18]), see for example FIG. |3] The 
most salient point is that this depth depends on kx and there- 
fore on E. Loosely speaking we can say that the further the 
energy is from £^ = the further it extends into the bulk. An 
upper limit for the momenta kx of the edge modes can be cal- 
culated from the condition that \di -\- < 1. Note that this 
condition also says that we must be inside the non-Abelian 
domain \J\ > 17^1/2 for the solution to be normalized. 
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FIG. 3. (Color online) (a) The function \f{yn) \ for different kx with Jz = 1, = —-7 and k = — .4. (b) A log plot of the same function 
fiVn) again with Jz = ^^J— —-7 and n — — A. (c) The penetration depth ^ as a function of kx for different values of J and fixed Jz and At. 
Penetration depth goes to infinity approximately when \di + d2\> 1 
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FIG. 4. A schematic of how 6 in the Majorana edge zero mode varies 
around an isolated domain of non-Abelian phase. In this model — 
tan~^(2Ai:/J). Inside the bulk, and at the comers, we indicate the 
phase that must be picked up as we move around that corner in the 
direction indicated by the arrows. 



V. FULLY OPEN BOUNDARY CONDITIONS 

If we surround a non-Abelian domain with an Abelian do- 
main we have no zero energy states if there are no vortices 
inside the non-Abelian domain. If we place an odd number of 
vortices inside the non-Abelian domain then we do have one 
zero energy edge mode even though an odd number of branch 
cuts intersect the domain wall. 

The key to understanding all this is that phases are also 
picked up when the wall direction is changed and that these 
phases all add up to tt, canceling the branch cut phase. A 
schematic of the phases picked up for the zero mode in a 
rectangular shaped system is shown in Figure |4] This pic- 
ture can be arrived at by analyzing each of the edges sepa- 
rately and assuming the appropriate plane wave momentum 



eigenstate ([16]) along each edge. The trigonometric identity 
tan~ la/6 + tan~ Ih/a = 7r/2 is the key to understanding 
why the total phase due to the corners is tt. At this time we 
have been unable to fully resolve the exact behavior at the cor- 
ners. However the numerically calculated example provided 
in Figure [4] shows that the phase changes at a corner happen 
abruptly and that the momentum eigenstates structure [16] is 
rapidly returned to as we move away from the comer. 

The chiral (non-zero energy) edge modes are also similar 
to that seen on the cylinder. In these cases, as for the zero 
momentum/energy modes, abrupt phase shifts are seen at the 
comers although in this case we cannot separate phase shifts 
due to momenta and those due to the corners. However it 
is worthwhile to note that if we use the value of momenta 
measured far from the corner in expression ( [Tt] ) we obtain 
the numerically calculated energy eigenvalue for the mode ex- 
actly. This measured value of momenta is however not exactly 
27rn/LTotai but slightly different magnitude. One could think 
of this as arising because the chiral mode sees a slightly dif- 
ferent perimeter Lrotai — ^l but we advise against taking this 
too literally. 

The picture above can be immediately applied to domains 
of Abelian phase inside a non-Abelian one. If there is no vor- 
tex inside this abelian domain then there is no branch cut and 
all modes are chiral but where the direction of the momenta 
for positive and negative energy modes is in the opposite sense 
to that on the outer edge. If an odd number of vortices exist 
inside the internal Abelian domain then we have an odd num- 
ber of branch cuts and a zero mode can exist. As suggested 
by Read and Green^ the zero mode due to a single vortex in 
the non-Abelian domain can be viewed as special case of this 
scenario where the domain edge has been reduced to a single 
plaquette. 



VI. CONCLUSION 

We have analysed edge mode structure of the Kitaev Hon- 
eycomb model using a Jordan-Wigner fermionization proce- 
dure. We see that the branch cuts are naturally defined for 
us with the single particle Hamiltonian ^ and the order pa- 
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(a) (b) (c) 



FIG. 5. (Color online) (a) The position space structure \uq\ of the Majorana edge zero-mode with a single vortex in the bulk (b) The phase 
dependence = ^ for the same Majorana zero mode (c) The phase dependence = 0-\-k-q for a. chiral edge mode where a vortex exists in the 
bulk. We can see here the combined effects of the branch cut, the phase jumps at corners and the almost constant momenta \k\ ^ 2n /LTotai- 
In these figures Jx — Jy — J — —0.7, Jz — 1, At = —.4 



rameter A. We then extended the notion of these branch cuts 
to account for edge effects between Abelian and non-Abelian 
domains. Although our general conclusions are in agreement 
with other methodolgies we feel there is an inherent simplic- 
ity to the above arguements that make them an important part 
of the overall story. 

For the specific model we have chosen we have been able 
to derive a simple recursive relation that exactly dictates the 
structure on edge between a vacuum and non-Abelian domain. 
A number of key features are present. Firstly the solutions 
are only normalized in the non-Abelian domain. Secondly 
we see a clear dependence on penetration depth on the mode 
momenta. We have also outlined how to apply the cylindrical 



solutions for the hard boundary to a fully open system. 

In future work we will attempt to analyse the edge mode 
momentum dependency further and to extend these results to 
softer boundaries. We will also attempt to identify enough 
properties to exactly formulate the mode structure at the cor- 
ners. 
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